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We discuss graviton oscillations based on the ghost free bi-gravity theory. We point out that 
this theory possesses a natural cosmological background solution which is very close to the case of 
general relativity. Furthermore, interesting parameter range of the graviton mass, which can be 
explored by the observations of gravitational waves, is not at all excluded by the constraint from 
the solar system tests. Therefore the graviton oscillation with possible inverse chirp signal would be 
an interesting scientific target of KAGRA, adv LIGO, adv Virgo and GEO. 



Introduction: Many works have been done for the 
detection possibility of modified propagation of gravita- 
tional waves due to finite graviton mass [HE]- However, 
adding mass to graviton was thought to be theoretically 
problematic due to the so-called Boulware-Deser (BD) 
ghost [3]. 

Recently, Hassan and Rosen proposed the first exam- 
ple of ghost-free bi-gravity models [4] , based on the fully 
nonlinear massive gravity theory in which the Boulware- 
Deser ghost is removed by construction [SHZ]. We con- 
sider two metrics expressed by 

•^•s = g^j.vdx^dx'^ , ds — (jf^i/dx^dx'^ . 

We introduce a ghost free action S = J d'^xC with 
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where Mq 



1/{8t:Gn)', Gn is the gravitational con- 
stant; YJf = y/gf^^gaiy', R and R are the Ricci scalars 
with respect to 5^^ and g^^, respectively; g and g are 
the determinants oi g^i, and g^jj, respectively; k is a con- 
stant which expresses the ratio between the two gravi- 
tational constants for 5^1, and g^^; Cn (n = 0, . . . , 4) are 
dimensionlcss constants, and Lm is the Lagrangian of the 
matter which interacts only with 5^^. By expressing the 
trace of F" as [F"] = tr(r") = Y^^^^Y;^^^ ■ • -Fa"""', we 
can write K,'s as 
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1, Vi = [Y], 
[Yf-3[Y][Y^]- 



V2 = [Y]' ~ [Y'] , 

-2[Y^], 
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The variation of the action with respect to g^^'^ and g'^'^ 
yields the field equations as 
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where T^^ is the energy momentum tensor of the ordinary 
matter, whereas 3^^ and B^^, come from the variations 
of the mass term. B^n and B^^, as well as T^^, satisfy 
conservation laws, which are explicitly given by 






0, V^T^'^0, V^B^-0, 



(1) 



where V and V are the covariant derivative operators 
with respect to g^j^y and g^^, respectively. 

The cosmological background: The background 
cosmology of this theory has been widely studied in 
Refs. [8-lOJ. but here our focus is on a particularly 
healthy branch. We assume that the two metrics can 
be written as 



ds^ =a^{-df 



dx^ 



dS^ ^d^{-c^df + dx^), 



where a, a, and c are functions of the time coordinate t. 
The Fricdmann equation for the physical metric reads 



3H^ = ^ 
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where we have introduced the Hubble parameter H = 
a I a}, the matter energy density pm including the dark 
energy, and the energy density due to the mass term 

py{i) = Mlm^ico + 3^ci + 6e'c2 + 6^03) , 



with ^ = a/a. The Friedmann equation for the hidden 
metric reads as follows 
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Writing down the first equation in Eq. M, we have 

3T{^)[daH ~{d./d)] =0, 

where r(^) = ci^ + 4c2^^ + 603^"^. This equation can be 
solved by imposing r(^) = or caH — (a/a) — 0, which 
implies the existence of two branches. In the following 
we will discuss the physical branch, defined by the latter 
condition, since the other branch is pathological^. Com- 
bining this condition with Eqs. (pi) and ([3]), we obtain an 
algebraic equation for ^ 
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If rrr 



>• Phi/Mq, the r.h.s of Eq. (4) should be very 
small. Denoting a value of ^ at which the right hand side 
vanishes by ^c, we focus on a cosmological background 
solution for which ^ asymptotes to ^c for pm — ^ 0. As we 
can absorb the constant part of pviO i^i^o the cosmolog- 
ical constant in pm, we also assume that pvi£,c) = 0. 

For this type of solution, we can expand ^ around ^c 
at low energies. Keeping only the linear order in ^ — i^c, 
Eq. Q becomes 
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where Fc = F(^c)- Substituting this relation into Eq. (pi), 
we recover the usual Friedmann equation as 
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with the effective gravitational constant given by 
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Mlil + nC). 



On using the definition of ^, the relation caH = a/ a 
implies S^ — {c — l)aH^. Substituting the differentiation 
of Eq. ^ into this relation, wc obtain 

c ~ 1 H ^7; — , 



^ The degrees of freedom of the theory reduce |10) . This will lead 
to a similar phenomenology observed in the original ghost-free 
single-metric massive gravity, which is characterized by the pres- 
ence of a scalar non-perturbative ghost 1111 . 



at low energies, where Pm is the matter pressure density. 
The above relation implies that the light cone of the hid- 
den metric automatically gets closer to the physical one 
as the matter energy density is diluted. 

Propagation of the gravitational waves: We now 
discuss the propagation of gravitational waves. We in- 
troduce tensor-type perturbations as gij = a'^{h+efj + 
hxe^j), and gij = a^ih+ef^ + hxS^j), with tr(e+e+) = 
1 ~ tr(£^e^), and tr(e+£^) = 0. The gravitational 
waves propagate at the speed of light for the physical 
sector, whereas at the speed c ~ 1 -I- 0{H^ /m?) for the 
hidden sector. However, the physical and hidden gravi- 
tons, because of the coupling through the mass term, will 
oscillate from one to the other. Keeping only the lead- 
ing effect of the deviation of c from unity, and neglecting 
the cosmic expansion effects, we write the propagation 
equations as [TD] 



h- /\h + m^Tcih ~h)=0, 
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where we have omitted the +/x index. For this set of 
equations, we write down the dispersion relation, assum- 
ing c — 1 ^ 1 but the magnitude of 

_ 2(27r/)^(g-l) 
is moderate, where we have defined 

-2_ (i + <?)rc™2 
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Then, for a given gravitational wave frequency /, two 
eigen wave numbers are given by 



fcL = (27r/)^ 
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and the corresponding eigen functions hi and ft.2 a-re re- 
lated to h and h as 



hi = cos.9gh -\- sinOg^/ K^c-h, 
/12 = — sin 9g h + cos 6g ^/k^c h 

with the mixing angle 
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We find that /i is the graviton mass of the second mode 
in the Minkowski limit (x — >■ 0). 

When we consider the propagation over a distance D, 
the phase shifts, due to the modified dispersion relation 
for their respective modes, are given by 
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FIG. 1: j<5<l>i,2| as a function of x for k^^ ~ 0.2 (dotted, 
black) , 1 (blue) and 100 (dashed, red) . Thick and thin curves 
represent |5$ij and J5<1>2|, respectively. 



Notice that this factor is symmetric under the replace- 
ment X — >■ l/x. In the limit a; — >■ 0, the first mode 
becomes massless. Although this mode also has non- 
trivial dispersion relation, its magnitude of modifica- 
tion tends to be suppressed. The factor fj,D\/c — 1 ~ 
\/3{l + kS,'^)VIqH D becomes 0(1) only after propagat- 
ing over a cosmological distance unless k^^ is extremely 
large, where VLq is the energy fraction of the dust matter 
at the present epoch. On the other hand, the remaining 
factor takes the maximum value 2~^/^ — (2 -(- 2k^^)~^/^ 
at a; = 1, which is also at most 0(1). In contrast to the 
first mode, the phase shift of the second mode can be 
significantly large when x is small or large. Here we plot 
5$i,2 in Fig. 1 for <2 = o.2, 1 and 100. 

Gravitational potential around a star in the 
Minkowski limit: In the above, we find that, unless 
K^^ is extremely large, a relatively small value of A^ to- 
gether with the excitation of the second mode is required 
for an observable magnitude of the phase shifts due to the 
non-trivial dispersion relation. Here we show that in the 
present bi-gravity models even with such a small value of 
A^ we can easily evade the solar system constraint from 
the precision measurement of gravity. 

In the low energy limit, it would be natural to assume 
the hierarchy, fc^ > /i^ > i/^. Since the Umit iJ -> 
is smooth, the _ff-dependent terms in the action appear 
as a positive power in H . Since such terms will not give 
any dominant contribution under the assumption of the 
above hierarchy, we set H — from the beginning here. 

Let us now consider static spherical symmetric per- 
turbations for both metrics induced by non-relativistic 
matter energy density /9,„, which is coupled only to the 
physical metric. We can write the respective perturbed 
metrics as 









without loss of generality. Here f is related to r by 
f = e^{r)r^ and Ti-ir) is another perturbation variable. 
We adopted the parametrization such that u vanishes in 
the case of general relativity. Now we write down the 
equations of motion and eliminate the variables on the 
hidden metric side, u, v and TZ. However, doing this is 
not so straightforward. In order to simplify the manipu- 
lation, we truncate the perturbation equations at second 
order and also neglect higher order terms in fj, appropri- 
ately. 

When we compute the terms second order in pertur- 
bation, we notice that there are terms enhanced by the 
factor l//i^. If we scrutinize these terms, some of them 
contain the factor 
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Our assumption here is that the energy scale of the bi- 
gravity theory itself is relatively high but the graviton 
mass /I is suppressed by a certain mechanism. Under this 
assumption, we pick up only the terms enhanced by the 
factor C / fi^ from the second order terms in the equations 
of motion. Then, after a little calculation, we obtain 
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where di is the differentiation with respect to the coor- 
dinates r(sin0cos(/>, sin0sin0, cos6') and A = 9^9% the 
standard three dimensional Laplacian operator, and we 
have defined C = C [1 -\- n(l) / {k(^) . At this level, the ex- 
pressions were recast into the form that does not assume 
spherical symmetry, where there is no ambiguity. 

Although we have truncated the equations at second 
order for simplicity, the higher order terms will not be 
suppressed once the second order terms become impor- 
tant. Nevertheless, such higher order terms will not 
change the following discussion as to the order of magni- 
tude estimate on the correction to the Newton's law. 

First we focus on Eq. (IT]) . Notice that non- vanishing u 
is the origin of the vDVZ discontinuity [T^] . This equa- 
tion tells us that the Vainshtein radius JIS] , within which 
the second term dominates the first term on the left hand 
side in Eq. ([7|, is given by 

ry = 0((Cr,A^)V3), 

where r^ is the gravitational radius of the star. From 
the above estimate, we find that the Vainshtein radius 
can be made arbitrarily large even with a large gravi- 
ton mass, if C is sufficiently large. Thus, the solar sys- 
tem can be easily contained within the Vainshtein ra- 
dius, where the second or even higher order terms on 
the left hand side of Eq. ^ dominate. Then, we have 
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Even if we require 



u to be smaller than 10 in the solar system r « 10 cm, 
A^ can be left arbitrarily small depending on the value 
ofC [14. 

Once u is suppressed on the scale of the solar sys- 
tem, Eq. (Is]) tells us that the equation for v does not 
largely deviate from the one in the Newtonian case: 
Aw = —Mq pm, and the gravitational constant is not 
different from the cosmological one. In Eq. (Is]) the terms 
second order in u are eliminated with the aid of Eq. ^ to 
make the effective gravitational source for v to be mani- 
fest. Solving the equation, we find that the correction to 
V is at most of 0{u). Notice that the mass term for v in 
Eq. (Is]) is absent. Therefore, v does not suffer from the 
Yukawa-type correction. 

Hence, one can conclude that the correction to the 

Newtonian potential is at most 0(^/ n^'^rgr/CX 
Namely, we can avoid the constraint from the test in 
the solar system, keeping the graviton mass sufficiently 
large. In the above we assumed that C is large. On the 
other hand, in constructing the cosmological background 
we have used the linear approximation for the deviation 
from the conformal equivalence between the two metrics, 
i.e. c— 1^1. If we further expand the background equa- 
tions in terms of c— 1, we find terms enhanced by the fac- 
tor C at second order. However, as long as CA^ < H~^, 
is satisfied, we can verify that the formulae for the back- 
ground metric remain approximately valid. In the early 
universe, where H is larger, the non-linear terms become 
necessarily important. However, the terms second order 
in ^ — ^c do not alter the effective Newton constant for 
the homogeneous background cosmology. 

The equations for u and v can be obtained similarly as 
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Once u is suppressed, i.e. if the Vainshtein mechanism 
is at work, we find v ~ v, which implies that metric 
perturbations on both sides are equally excited by the 
matter fields. 

Graviton oscillations and inverse chirp sig- 
nal: Here we begin with discussing the generation of 
gravitational waves. We found that the metric excita- 
tions are almost conformal within the Vainshtein radius 
of a star. If we consider the junction between the near- 
zone metric perturbation with the far-zone metric de- 
scribed as gravitational waves, both h and h are excited 
exactly as in the case of general relativity. This implies 
that both eigen modes hi and /12 are excited unless x — 0. 
(Recall that h2 oc h — h when a; = 0.) 

One may suspect that the linear approximation to the 
gravitational wave perturbation equations ^ is not valid 
within the Vainshtein radius. However, the effective en- 
ergy momentum tensor coming from the variation of the 
mass term, which gives corrections to the case of general 
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FIG. 2: _Bi,2 as a function of x for k^^ = 0.2 (dotted, black), 
l(blue) and 100(dashed, red). Thick and thin curves represent 
Bi and B2, respectively. 



relativity, is largely enhanced only for the terms purely 
composed of w (or equivalently u), which behave as clouds 
around localized matter sources. Namely, it just con- 
tributes as the source of gravitational waves but does 
not change the wave propagation. The other corrections 
are suppressed as long as the amplitude of the deviation 
of the metric from the case of general relativity remains 
small. 

Next, we analyze the gravitational waveform from in- 
spirals of NS-NS binaries at a distance. For the current 
bi-gravity model, our detector signal becomes a linear 
combination of two components, whose relative ampli- 
tudes are determined by the mixing angle 9g. For sim- 
plicity, we here neglect the time dependence of 9g as well 
as all the cosmological effects. Using the stationary phase 
approximation and flux conservation, the observed signal 
is given in Fourier space as 



h{f) = A(/)e'*(-^) iBie'^^^'-f^+B^e 
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where the amplitude A{f) (after angular average), B12 
and the phase function $(/, 5) (truncated at 1.5PN or- 
der) are given by 



A{f) = 
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30 D 

El = cos 9g (cos 6g + a/k^c sin 9g ) , 
B2 = sin 9g (sin 9g — a/k^c cos 9g ) , 
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with y = TrMf, the chirp mass A4 = 
(77111712)^'^ /{mi + 7712)^'^ and the reduced mass ra- 
tio r] = mim2/{mi + m2)'^. The first and second 
terms in Eq. (J9]) show the contributions of hi and /12, 
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FIG. 3: The arrival time as functions of the frequency / 
for respective modes for IAMq + IAMq binary inspiral with 
<^ = 100, D = SOOMpc H = 67.3km s-^Mpc'S tlo = 0.315 
and Ap = O.OOlpc. The blue solid curve is for the first mode, 
while the dashed red one for the second mode. 



respectively. Here we plot B12 in Fig. 2 for k^^ — 0.2, 1 
and 100. 

For a; ^ 1, the excitation of the second mode /12 is sup- 
pressed. Furthermore, 54>i is suppressed in this regime. 
Therefore, the propagation of gravitational waves is sim- 
ilar to the case of general relativity. For x ^ 1, both hi 
and /12 are equally excited. However, since the gravita- 
tional wave detector can detect the perturbation of the 
physical metric only, we can observe only hi. Therefore, 
the frequencies at which both modes can be observed are 
limited to a; w 1. This is the meaning of Fig. 2. 

When both modes are observable, graviton oscillations 
due to the interference between two modes can be de- 
tected beyond the distance scale where 5<^2 ~ S'^i be- 
comes 0(1). The difference of the phases 5(t)i — (5(/)2 
is minimum at x = 1 as shown in Fig. 1, which is 
evaluated as 5(f)i — 5(t>2\x=i — ^/MTqHD. Therefore, 
one may think that the effect is really small as long as 
D < H~^ . However, the average density of the uni- 
verse is much lower than the average density in galax- 
ies, where binaries are embedded. Therefore gravita- 
tional waves experience much lower value of cc, typically 
X « 10~®, during the propagation. Roughly speaking, 
<5$2 - (5$i ~ v^3(l + K\^)no/2xHD for x < 1. Hence, 
the effect can be largely enhanced. 

Once d{6^2 — 5^i)/df becomes sufficiently large, the 
arrival times of two modes are different. Then, we 
may observe two chirp signals. Using the stationary 
phase approximation, the relation between the arrival 
time of the wave and the frequency is determined by 
t — d{^ + S^i)/d{2Trf). As an illustrative purpose, 

^ In fact, another possibility of suppressing u in Eq. ¥n is to con- 
sider a really large value such as /i~^ < O.lmni, In this case 
the massive graviton will be almost irrelevant to describe the 



in Fig. 3 we show the shifts of the arrival time com- 
pared with the case of general relativity for k£,'^ = 100, 
D = aOOMpc H = 67.3km s-^Mpc"!, f^o = 0.315 and 
fi = (0.001pc)-\ for which a; « 4 x 10"** at lOOHz. One 
can see that the relation between the arrival time and 
the frequency is reversed for the second mode, i.e. in- 
verse chirp signal may occur. 

Since a large graviton mass can be consistent with the 
solar system test in our current modeP, there is a possi- 
bility that we may detect the graviton oscillations or even 
the inverse chirp signal by the next generation gravita- 
tional wave detectors. In the present model, measurable 
effects are expected only when k^^ is large. This require- 
ment may cause some conflict with observations but at 
flrst analysis there seems to be no severe constraint. We 
think this model gives the first existence proof of models 
in which a measurable deviation from general relativity 
in the gravitational wave propagation can be expected. 
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